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Abstract 

We review some recent results on asymptotic properties of polynomials of 
large degree, of general holomorphic sections of high powers of positive line 
bundles over Kahler manifolds, and of Laplace eigenfunctions of large eigen- 
value on compact Riemannian manifolds. We describe statistical patterns in 
the zeros, critical points and L p norms of random polynomials and holomor- 
phic sections, and the influence of the Newton polytope on these patterns. 
For eigenfunctions, we discuss L p norms and mass concentration of individual 
eigenfunctions and their relation to dynamics of the geodesic flow. 
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1. Introduction 

In many measures of 'complexity', eigenfunctions V~K<p\ = \tp\ of first order 
elliptic operators behave like polynomials p(x) = Y^,\ a \<x c a x ° of degree N ~ A [6]. 
The basic example we have in mind is the Laplacian A on a compact Riemannian 
manifold (M, g), but the same is true of Schroedinger operators. The comparison 
is more than an analogy, since polynomials of degree N are eigenfunctions of a first 
order elliptic system. 

The comparison between eigenfunctions and polynomials is an essentially local 
one, most accurate on small balls B(xo, j). Globally, eigenfunctions reflect the 
dynamics of the geodesic flow G* on the unit (co-)tangent bundle S*M . This is one 
of the principal themes of quantum chaos. 
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In this article, we review some recent results on the asymptotics of polynomials 
and cigcnfunctions, concentrating on our work in collaboration with P. Blchcr, A. 
Hassell, B. Shiffman, C. Sogge, J. Toth and M. Zworski. A unifying feature is the 
asymptotic properties of reproducing kernels, namely Szcgo kernels ILn(z,w) in the 
case of polynomials, and spectral projections E\(x, y) for intervals [A, A + 1] in the 
case of eigenfunctions of \fK. For other recent expository articles, see [9, 26]. 

2. Polynomials 

There are several sources of interest in random polynomials. One is the desire 
to understand typical properties of real and complex algebraic varieties, and how 
they depend on the coefficients of the defining equations. Another is their use as a 
model for the local behavior of more general eigenfunctions. A third is that they 
may be viewed as the eigenvectors of random matrices. Just as random matrices 
model the spectra of 'quantum chaotic' systems, so random polynomials model their 
eigenfunctions. 

2.1. SU(m + 1) polynomials on CP m and holomorphic sec- 
tions 

Complex polynomials of degree < p in m variables form the vector space 
V™:={f(z u ...,z m )= Yl c a z^---z%T, c a eC}. 

a£N m :\a\<p 

To put a probability measure on V™ is to regard the coefficients c a as random 
variables. The simplest measures are Gaussian measures corresponding to inner 
products on "P™. By homogenizing / to F(zo, z\, ■ ■ ■ ,z m ) of degree p, we may 
identify V™ with the space H°(CP m , 0(p)) of holomorphic sections of thepth power 
of the hyperplane bundle. It carries the standard SU(m + l)-invariant Fubini- 
Study inner product (F\,F2)fs = J S 2 m +i FiF^da , where da is Haar measure on 
the (2m + l)-sphere S 2m+1 . An orthonormal basis of ff°(CP m , 0(p)) is given by 
{ m an s }. The corresponding SU(m + l)-invariant Gaussian measure 7,5 is defined 

dls{s) = i eHA|2dA ' s= £ Xa w^- 

\a\<p 

Thus, the coefficients A a are independent complex Gaussian random variables with 
mean zero and variance one. 

More generally, we can define Gaussian ensembles of holomorphic sections 
H°(M, L N ) of powers of a positive line bundle over any Kahler manifold (M,oj). 
Endowing L with the unique hermitian metric h of curvature form oj, we induce 
an inner product (,) on H°(M,L N ) and a Gaussian measure jn. We denote the 
unit sphere in H°(M, L N ) relative to this inner product by SH°(M, L N ). The Haar 
measure on SH°(M, L N ) will be denoted ^n- It is closely related to the Gaussian 
measure. 
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2.2. 



Zeros 



The problems we discuss in this section involve the geometry of zeros of sections 
s £ H°(M, L N ) of general positive line bundles. There is a similar story for critical 
points. 

• Problem 1 How are the simultaneous zeros Z s = {z : Si(z) = ■ ■ ■ = Sk(z) = 
0} of a k-tuple s = (si, . . . , Sfc) of typical holomorphic sections distributed? 

• Problem 2 How are the zeros correlated? When k = m, the simultaneous 
zeros form a discrete set. Do zeros repel each other like charged particles? Or 
behave independently like particles of an ideal gas? Or attract like gravitating 
particles? 

By the distribution of zeros we mean either the current of integration over Z s or 
more simply the Riemannian (2m— 2fc)-volume measure (|Z S |, ip) = J z (pdVol^m-^k ■ 



where \Z s \ n denotes the product of the measures \Z S \ on the punctured product 
M„ = {(z 1 , . . . , z n ) e M x • • • x M : z p ^ z q for p ^ q} and where dz denotes the 
product volume form on M n . 

The answer to Problem 1 is that zeros almost surely become uniformly dis- 
tributed relative to the curvature u of the line bundle [18]. Curvature causes sec- 
tions to oscillate more rapidly and hence to vanish more often. More precisely, we 
consider the space S = IIiv=i SH°(M 7 L N ) of random sequences, equipped with 
the product measure measure \i = IIn^iMw- An element in S will be denoted 
s = {sw}. Then, jjZ s — > ui, as N — > oo for almost every s. 

The answer to Problem 2 is more subtle: it depends on the dimension. We 
assume fc = m so that almost surely the simultaneous zeros of the fc-tuple of sections 
form a discrete set. We find that these zeros behave almost independently if they 
are of distance > apart for D ^> 1. So they only interact on distance scales 

of size ^t=. Since also the density of zeros in a unit ball Bi(zq) around zq grows 

like N m , we rescale the zeros in the 1/%/iV-ball B l ^^(zo) by a factor of V~N to 
get configurations of zeros with a constant density as N — > oo. We thus rescale the 
correlation functions and take the scaling limits 



In [1] , we proved that the scaling limits of these correlation functions were univer- 
sal, i.e. independent of M,L,uj,h. They depend only on the dimension m of the 
manifold and the codimension k of the zero set. 

In [2], we found explicit formulae for these universal scaling limits. In the case 
n = 2, K^j!. m (z 1 , z 2 ), depends only on the distance between the points z 1 , z 2 , since 
it is universal and hence invariant under rigid motions. Hence it may be written as: 





z n )= Jim K? k (zo)- n K» k (zo + 4z 




(2.1) 



) = Kkmilz 1 - Z 2 \) . 



(2.2) 



We refer to [1] for details. 
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Theorem 1 [2] The pair correlation functions of zeros when k = m are given by 
[ m±l r 4 - 2m + 0(r s - 2m ) , as r^O 

Kmm(r) = < (2.3) 

{ l + 0(e- Cr ), (C>0) asr^oo. 

When m = l,K mm (r) — > as r — > and one has "zero repulsion." When 
m = 2, K mm (r) — > 3/4 as r — ► and one has a kind of neutrality. With m > 3, 
n m m(r) / oo as r — > and there is some kind of attraction between zeros. More 
precisely, in dimensions greater than 2, one is more likely to find a zero at a small 
distance r from another zero than at a small distance r from a given point; i.e., 
zeros tend to clump together in high dimensions. 

One can understand this dimensional dependence heuristically in terms of 
the geometry of the discriminant varieties V 1 ^ C H°(M 7 L N ) m of systems S — 
(si, . . . ,s m ) of m sections with a 'double zero'. The 'separation number' sep(F) 
of a system is the minimal distance between a pair of its zeros. Since the nearest 
element of V 7 ^ to F is likely to have a simple double zero, one expects: sep(F) <~ 
^dist{F,V^). Now,the degree of V™ is approximately N m . Hence, the tube (V%) e 
of radius e contains a volume ~ e 2 N m . When e <~ jV~ m / 2 , the tube should cover 
PH°(M, L N ). Hence, any section should have a pair of zeros whose separation is 
<~ N~ m / 4 apart. It is clear that this separation is larger than, equal to or less than 
iV -1 / 2 accordingly as m = 1, m = 2, m > 3. 

2.3. Bergman-Szego kernels 

A key object in the proof of these results is the Bergman-Szego kernel Hn(x, y), 
i.e. the kernel of the orthogonal projection onto H°(M,L N ) with respect to the 
Kahlcr form u>. For instance, the expected distribution of zeros is given by Ejv(Z/) = 
^3!<991ogn7v(z, z) + u). Of even greater use is the joint probability distribu- 
tion (JPD) Dn(x 1 , . . . ,x n ;^ 1 , . . . ,£ n ; z 1 , . . . ,z n ) of the random variables x^{s) — 
s(zi), &(s) — Vs(z J ), which may be expressed in terms of II jv and its deriva- 
tives. In turn, the correlation functions may be expressed in terms of the JPD by 

K N (z\...,z n ) = jD N (0^,z)UU(UWH t 1 ]-. 

The scaling asymptotics of the correlation functions then reduce to scaling 

asymptotics of the Bergman-Szego kernel: In normal coordinates {zj} at Jo € M 

and in a 'preferred' local frame for L, we have [1]: 

_n^(p + — ,-,p + — ,-) 

_ e i { e- v ) + u-v-U\u\ 2 +\v\-) j! + &i(U} v)N -\ + . . . - _ 

To be precise, IT^v is the natural lift of the kernel as an equivariant kernel on the 
boundary dD* of the unit (co-) disc bundle of L* . Note that e ^-f)+u-v-^(\u\ 2 +\v\ 2 ) 
is the Bergman-Szego kernel of the Heisenberg group. These asymptotics use the 
Boutet de Monvel -Sjostrand parametrix for the Bergman-Szego kernel [4], as ap- 
plied in [29] to the Fourier coefficients of the kernel on powers of positive line 
bundles. 
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2.4. Polynomials with fixed Newton polytope 

The well-known Bernstein-Kouchnircnko theorem states that the number of 
simultaneous zeros of (a generic family of) m polynomials with Newton polytope 
P equals m\Vol(P). Recall that the Newton polytope Pf of a polynomial is the 
convex hull of its support Sf = {a € Z m : c a ^ 0}. Using the homogenization map 
/ — > F, the space of polynomials / whose Newton polytope Pf contained in P may 
be identified with a subspace 

H°(CP m , 0(p), P) = {Fe H°(CP rn , 0(p)) -PfCP} (2.4) 

of H (CP m ,O{p)). 

The problem we address in this section is: 

• Problem 3 How does the Newton polytope influence on the distribution of 
zeros of polynomials? 

Again, one could ask the same question about L 2 mass, critical points and so 
on and obtain a similar story. In [19] we explore this influence in a statistical and 
asymptotic sense. The main theme is that for each property of polynomials under 
study, P gives rise to classically allowed regions where the behavior is the same as 
if no condition were placed on the polynomials, and classically forbidden regions 
where the behavior is exotic. 

Let us define these terms. If P C is a convex integral polytope, then the 
classically allowed region for polynomials in iJ°(CP m , 0(p), P) is the set 

A P := Ms 1 C C* m 

(where P° denotes the interior of P), and the classically forbidden region is its 
complement C* m \Ap. Here, ^(z) = (yqq^p, • • • , i+fjjp) is the moment map of 
CP" 1 . 

The result alluded to above is statistical. Since we view the polytope P of 
degree p as placing a condition on the Gaussian ensemble of SU(P) polynomials 
of degree p, we endow H°(CP m 7 0(p), P) with the conditional probability measure 
ls\p- 

dls\p{s) = ^pe-W d\, S =^A Q ^-, (2.5) 

a£P 11 1 

where the coefficients A Q are again independent complex Gaussian random variables 
with mean zero and variance one. 

Our simplest result concerns the the expected density Eip(Z/ 1 f m ) of the 

simultaneous zeros of {fi, ■ ■ ■ , f m ) chosen independently from H°(CP m , 0{p) 1 P). 
It is the measure on C* m given by 

E|p(%,...,/J(t0 

= j dj plP (h) ■ ■ ■ j d lplP (f m ) [#{z e U : fi(z) = ■ ■ ■ = f m (z) - 0}] , (2.6) 
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for U C C* m , where the integrals are over H°(CP m , 0{p), P). We will determine 
the asymptotics of the expected density as the polytope is dilated P — > NP, JVeN. 

Theorem 2 [19] Suppose that P is a simple polytope in R m . Then, as P is dilated 
toNP, 

I [ Wpg on Ap 



on C* m \.4 



p 



in the distribution sense; i.e., for any open U C C* m , we have 

1 -E| NP (#{z e U : fi{z) = ■ ■ ■ = f m (z) = 0}) -> m!Vol Cpm ([/n A P ) 



(Nsy 

There are also results for k < m polynomials. The distribution of zeros is Wp S 
in Ap as if there were no constraint, while there is an exotic distribution in C*" 1 \ 
Ap which depends on the exponentially decaying asymptotics of the conditional 
Bergman- Szego kernel 

Ez a w a 
r, — n n n — j 

i.e. the orthogonal projection onto the subspace (2.4). It is obtained by sifting 
out terms in the (elementary) Szego projector of H°(CP m , 0(pN)) using the poly- 
tope character XNp{e lip ) — J2 aeNP e 1 ^'^ . To obtain asymptotics in the forbidden 
region, we write XNp{e tip ) = J Mp Tl^ p (e t ' p w,w)dV(w), where n Mp is Bergman- 
Szego kernel of the toric variety Mp associated to P. We then make an explicit 
construction of H^ F as a complex oscillatory integral. An alternative is to express 
Xnp as a Todd derivative of an exponential integral over P (following works of 
Khovanskii-Pukhlikov, Brion-Vcrgnc and Guillemin). We thus obtain a complex 
oscillatory integral formula for U\ NP (z, w). To obtain asymptotics in the forbid- 
den region we carefully deform the contour into the complex and apply a complex 
stationary phase method. 

Although wc only discuss expected behavior of zeros here, the distribution of 
zeros is self-averaging: i.e., almost all polynomials exhibit the expected behavior in 
an asymptotic sense. We also expect similar results for critical points. 



3. Eigenfunctions 

We now turn to the eigenvalue problem ^g(p^ — X^tpv, (tp v ,Lp v i^ — o~w' on a 
compact Riemannian manifold (M,g). We denote the A-eigenspace by V\. The role 
of the Szego kernel is now played by the kernel E\(x,y) = J2x Kx^vi^fviv) °f 
the spectral projections. 



3.1. L p bounds 



Asymptotics of Polynomials and Eigenfunctions 



739 



Our first concern is with LP norms of L 2 -normalized eigenfunctions. We mea- 
sure the growth rate of LP norms by LP (A, g) = sup^gy^.n^n 2=1 ||</?||£ P . By the local 
Weyl law, E x (x,x) = Ea„<a \<pA x )\ 2 = ( 27 0~" I\ 6 \<x d t + 0(A n ~ 1 ), it follows that 

L°°(X,g) = 0(A Ii 2 i ) on any compact Riemannian manifold. This bound, which is 
based entirely on a local analysis, is sharp in the case of the standard round sphere, 
S n or on any rotationally invariant metric on S 2 , but is far off in the case of flat 
tori. This motivates: 

• Problem 5 For which (M,g) is this estimate sharp? Which (M,g) are ex- 
tremal for growth rates of ||<^a||p> both maximal and minimal? What if M 
has a boundary? What is the expected LP norm of a 'random' L 2 -normalized 
polynomial or eigenfunction? 

In [20] , we give a necessary condition for maximal eigenfunction growth: there 
must exist a point x € M for which the set C x = {£ e S*M : 3T : exp x T£ = x} of 
directions of geodesic loops at x has positive surface measure. 



Theorem 3 [20] Lf C x has measure in S*M for every x £ M then 

(3.1) 



L?(\, 9 ) = o(X^), P >*W *M = {2LiV 

[ 2 '2 p/> 



2 < p < M . 

— r — n—l 



The L p -bounds 0(\ s{:p) ) were proved by Sogge to hold for all (M,g). 
We further prove: 

Theorem 4 [20] (see also [17]) Suppose that (M,g) is: 

• Real analytic and that L°°{\g) = ^(A^- 1 )/ 2 ). Then (M,g) is a Y £ m - 
manifold, i. e. 3m such that all geodesies issuing from the point m return to 
m at time I. Ln particular, if dim M — 2, then M is topologically a 2-sphere 
S 2 or a real projective plane M.P 2 . 

• Generic. Then L°°(A, 5 ) = o(A(™- 1 )/ 2 ). 

Here, fi(A~) means O(At) but not The generic result holds because 

C x has measure in S*M for all x e M for a residual set of metrics. 

In the case of random polynomials, or random combinations of eigenfunctions 
in short spectral intervals, the almost sure growth of L°° norms is 0(\/log N) while 
the LP norms for p < 00 are bounded. This was proved by J. Vanderkam [24] for 
S m , Nonnenmacher-Voros [14] for elliptic curves and Shiffman-Zelditch (to appear) 
for the general case using Levy concentration of measure estimates. 



3.2. Integrable case 

Results on minimal growth have been obtained by J. A. Toth and the author 
in the quantum completely integrable case, where y/A = P\ commutes with n—l 
first order pseudodifferential operators P2, ■ . ■ , P n € 5 ,1 (M) (n = dim M) satisfying 
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[Pi, Pj] = and whose symbols define a moment map V := (pi, . . . ,p n ) satisfying 
dpi A dp2 A • • • A dp n 7^ on a dense open set C T*M — 0. Since {pi 7 pj} = 0, the 
functions pi,...,p n generate a homogeneous Hamiltonian Reaction whose orbits 
foliate T*M — 0. We refer to this foliation as the Liouville foliation. 

We consider the LP norms of the L 2 -normalized joint eigcnfunctions Pjfx = 
Xjipx. The spectrum of A often has bounded multiplicity, so the behaviour of joint 
eigenfunctions has implications for all eigenfunctions. 

Theorem 5 [22, 23] Suppose that the Laplacian A g of(M,g) is quantum completely 
integrable and that the joint eigenfunctions have uniformly bounded L°° norms. 
Then (M, g) is a flat torus. 

This is a kind of quantum analogue of the 'Hopf conjecture' (proved by Burago- 
Ivanov) that metrics on tori without conjugate points are flat. In [23], a quantitative 
improvement is given under a further non-degeneracy assumption. Unless (M, g) is 
a flat torus, the Liouville foliation must possess a singular leaf of dimension < n. 
Let £ denote the minimum dimension of the leaves. We then construct a sequence 
of eigenfunctions satisfying: 

n-t (n-l)(p-2) 

Ikfclkoc >C(e)A fc 4 , \\<fik\\LP>C(e)X k 4p , (2 < p) 

for any e > 0. It is easy to construct examples were £ — n—1, but it seems plausible 
that in 'many' cases £ = 1. To investigate this, one would study the boundary faces 
of the image V(T* M — 0) of T * M — under a homogeneous moment map. For a 
related study in the case of torus actions, see Lerman-Shirokova [12]. 

3.3. Quantum ergodicity 

Quantum ergodicity is concerned with the sums (A G \I ,0 (M)): 

S P (X)= Yl \(Aip v ,<p v )-u{A)\>, w(A)= / OAdn. (3.2) 

In work of A.I. Schnirclman [11], Colin de Verdiere and the author [27], it is 
shown that S P (X) = o(N(X)) if G* is ergodic. In the author's view [27], this is best 
viewed as a convexity theorem. We mention briefly some new results. 

In work of Gerard-Leichtnam [7] and Zelditch-Zworski [30], the ergodicity re- 
sult was extended to domains dfl with piecewise smooth boundary and ergodic 
billiard flow. Since the billiard map on B*dQ is ergodic whenever the billiard flow 
is, suitable boundary values of ergodic eigenfunctions (e.g. <pk\dn in the Neumann 
case or d v <fk\dn m the Dirichlct case) should also have the ergodic property. This 
was conjectured by S. Ozawa in 1993. A proof is given in our work with A. Has- 
sell [8] for convex piecewise smooth domains with ergodic billiards (in the case of 
domains with Lipschitz normal and with Dirichlet boundary conditions, this had 
earlier been proved in [7] by a different method). 

Little is known about the rate of decay. For sufficiently chaotic systems 
(satisfying the central limit theorem) one can get the tiny improvement S P {X) = 
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0(A(A)/(logA) p )) [28]. The asymptotics S 2 {\) ~ B(A)X have recently been ob- 
tained by Luo-Sarnak [13] for Hecke eigenfunctions of the modular group, exploiting 
the connections with L-functions. These asymptotics (though not the coefficient) 
are predicted by the random polynomial model. Other strong bounds in the arith- 
metic case were obtained by Kurlbcrg-Rudnick for eigenfunctions of certain quan- 
tized torus automorphisms [10]. Bourgain-Lindenstrauss [3] and Wolpert [25] have 
developed the 'non-scarring' result of [16] to give entropy estimates of possible 
quantum limit measures in arithmetic cases. 
A natural problem is the converse: 

• Problem 6 What can be said of the dynamics if S P (X) = o(N(\))l Docs 
quantum ergodicity imply classical ergodicity? 

It is known that classical ergodicity is equivalent to this bound plus estimates 
on off-diagonal terms [21]. The existence of KAM quasimodes (due to Lazutkin 
[11], Colin de Verdiere [5], and Popov [15]) makes it very likely that KAM systems 
are not quantum ergodic, nor are (M, g) which have stable elliptic orbits. 

A further problem which may be accessible is: 

• Problem 7 How are the nodal sets {</?„ = 0} distributed in the limit v — > oo? 

In [14] (for elliptic curves) and [19] (general Kahler manifolds) it is proved that 
the complex zeros of quantum ergodic eigenfunctions become uniformly distributed 
relative to the volume form. Can one prove an analogue for the real zeros? 
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